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Reproducible scanning tunneling microscope (STM) spectra of BiiSr-zCaCu-zOs+s consistently ex- 
hibit asymmetric tunneling characteristics, with the higher peak conductance corresponding to a 
negatively biased sample. We consider various possible sources of this asymmetry that are not intrin- 
sic to the superconducting state, including energy dependence of the normal state densities of states 
of sample and/or tip, existence of bandwidth cutoffs, unequal work functions of tip and sample, and 
energy-dependent transmission probability. It is concluded that none of these effects can explain 
the sign and temperature dependence of the observed asymmetry. This indicates that the observed 
asymmetry reflects an intrinsic property of the superconducting state: an energy-dependent super- 
conducting gap function with non-zero slope at the Fermi energy. It is pointed out that such a sloped 
gap function will also give rise to a thermoelectric effect in STM experiments, resulting in a positive 
thermopower. We discuss the feasibility of observing this thermoelectric effect with an STM and 
conclude that it is easily observable. Again, contributions to this thermoelectric effect may also arise 
from energy dependence of normal state densities of states and from energy-dependent transmission 
probability. However, because each of these features manifests itself differently in the thermoelec- 
tric effect and in the tunneling characteristics, an analysis of thermoelectric currents and voltages 
together with the tunneling spectra as function of temperature and tip-sample distance would allow 
for accurate determination of the slope of the gap function. It is suggested that it would be very 
worthwhile to perform these experiments, because the slope of the gap function reflects a fundamen- 
tal property of the superconducting state. In particular, the theory of hole superconductivity has 
predicted the existence of such a slope, of universal sign, in all superconductors. It is furthermore 
argued that recent experimental results on vortex lattice imaging provide further strong evidence 
for the existence of the gap slope discussed here. 



I. INTRODUCTION 

Tunneling spectroscopy is a powerful experimental tool 
that has the potential to provide fundamental informa- 
tion on the microscopic physics of high temperature su- 
perconductivity. Unfortunately, due to the difficulty in 
preparation of well characterized sample surfaces and 
tunnel junctions, often the spectra obtained have shown 
varying features that depend on the particular sample 
and/or the particular tunneling technique used. Conse- 
quently it has been difficult to determine which of the 
features found in the large number of experimental stud- 
ies that have been performed to date actually reflect in- 
trinsic properties of high temperature superconductors. 

It is hence notable that several recent scanning tun- 
neling microscopy (STM) studies of single crystals of 
BiiSriCaCu20%+b cleaved in ultrahigh vacuum have 
shown highly reproducible features. [PH§ These exper- 
iments arc performed under carefully controlled condi- 
tions, and tunneling spectra are obtained at a large num- 
ber of different positions on the sample and for a range 
of distances between tip and sample. The spectra do not 
show time dependence and generally exhibit either super- 
conducting or semiconducting features, which has been 



suggested to depend on the oxygen stochiometry of the 
topmost Cu-0 layer at the position of the tip Q . 

The spectra which exhibit superconducting features 
show structure in the tunneling conductance dl/dV that 
resembles the structure in the tunneling density of states 
of conventional BCS superconductors predicted by BCS 
theory: a decrease in the conductance for low voltages, 
where a gap-like structure develops, and a large enhance- 
ment of the conductance at the edges of this gap-like 
structure. A notable difference with the predictions of 
conventional BCS theory is that the spectrum is highly 
asymmetric: the peaks in dl/dV are symmetrically lo- 
cated around zero voltage but their height is substan- 
tially different, with the largest peak corresponding to 
a negatively biased sample. This asymmetry is found 
to persist to low temperatures. Another feature of the 
spectra, which is sometimes observed only for negative 
sample voltage JjJ , and sometimes for both polarities || , 
is a dip in the conductance at voltages about twice the 
position of the conductance peaks. 

In this paper we focus on the tunneling asymmetry in 
the peaks of dl/dV and consider various possible origins 
for it. Firstly, it is clear that the tunneling geometry 
is asymmetric itself; an asymmetry may be expected to 
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arise from different work functions of the tip and the 
sample. In fact, tunneling asymmetries are often seen in 
conventional tunnel junctions and have been explained 
assuming differences in work functions of the electrodes 
P]. Secondly, as we will show, a tunneling asymmetry 
also arises from the fact that the tunneling probability 
is energy- dependent, even in the absence of differences 
in work functions. Finally, a tunneling asymmetry could 
also arise from an energy-dependent normal state density 
of states of either tip or sample, or from the existence of 
a bandwidth cutoff in cither tip or sample. To focus 
the discussion we assume that the superconducting state 
in the sample is described by conventional BCS theory 
with an isotropic gap, so that the density of states in the 
superconducting state is 



N(E(e)) = g s (e) 



E 



(1) 



with g s (e) the density of states of the sample in the nor- 
mal state, and the quasiparticle energy E given in terms 
of the band energy e by 



(2) 



with A the energy gap and /i the chemical potential. We 
will also allow for the possibility of broadening of the 
density of states arising from lifetime or other effects || , 
characterized by an inverse lifetime T. 

We will show that none of the effects listed above is 
able to describe even qualitatively an asymmetry of the 
sign and temperature dependence seen experimentally. 
Briefly, asymmetries originating in asymmetric normal 
state densities of states or bandwidth cutoffs tend to be- 
come smaller as the temperature is lowered, contrary to 
experimental observations. Asymmetry originating in en- 
ergy dependence of the transmission probability persists 
down to low temperatures but is of opposite sign to that 
observed experimentally. Asymmetry originating in dif- 
ferences in work functions can give the correct sign of the 
observed asymmetry at high temperature but the sign 
switches as the temperature is lowered, in contradiction 
with experiment. 

We are thus led to the conclusion that the observed 
asymmetry originates in an intrinsic feature of the su- 
perconducting state. Such asymmetry will arise when 
the superconducting gap function A is itself energy de- 
pendent Q , so that the quasiparticle energy Eq. (2) is 



E{e) = y/(e-tf + A(e)*. 



(3) 



The measured asymmetry reflects the slope of the gap 
function at the Fermi energy, A'(/i). Hence we argue 
that the experimental results are direct evidence for the 
existence of a slope in the superconducting gap function. 

Another experimental consequence of a sloped gap 
function is that it should manifest itself in a thermo- 
electric effect in tunneling experiments. It gives rise 
to a thermoelectric power of the tunnel junction given by 



Q — A(^)A'(/i)/eT for small temperature gradients. We 
argue that such a thermoelectric effect should be easily 
observable with an STM and that it should be experi- 
mentally searched for. Certainly, energy dependence of 
the normal state densities of states and energy-dependent 
transmission could also contribute to a thermoelectric ef- 
fect. However we will show that analysis of the ther- 
moelectric currents and voltages in conjunction with the 
tunneling characteristics dl/dV as function of tempera- 
ture and tip-sample distance should allow to clearly dis- 
tinguish between the different sources of a thermoelec- 
tric effect, and in particular to accurately determine the 
slope of the gap function. The sign of the observed tun- 
neling asymmetry implies that the sign of the measured 
thermoelectric power will be positive (at least for small 
tip-sample distance), opposite to what one would pre- 
dict from energy dependence of the electron transmission 
probability alone. In fact, thermoelectric experiments 
with STM with normal metal electrodes have recently 
been reported with the sign of the measured the- 
mopower always negative. 

The slope of the gap function at the Fermi energy is a 
fundamental property of the superconducting state, as it 
gives information on the underlying microscopic interac- 
tions that give rise to superconductivity in the system. 
For example, for an attractive Hubbard model [[| or for 
an electron-phonon model within a local approximation 
(Holstcin model) 10 the superconducting gap function 
has zero slope. The detection of a gap slope would re- 
veal the existence of an intrinsic electron-hole asymmetry 
of the superconductor and shed light on the theoretical 
understanding of superconductivity. In particular, the 
theory of hole superconductivity has predicted the 
existence of such a slope, of universal sign, in all super- 
conductors, with magnitude that scales with the criti- 
cal temperature. However, a sloped superconducting gap 
function could also arise in other theoretical frameworks. 

It is currently generally believed that high T c super- 
conductors have an order parameter with complete or 
dominant d-wave symmetry |12|. The d-wave symmetry 
relates to the variation of the gap on the Fermi surface, 
while the slope of the gap function discussed here relates 
to variation of the gap in direction perpendicular to the 
Fermi surface. Here we assume a gap that is constant on 
the Fermi surface for simplicity; however a simple gen- 
eralization of our analysis could also be used to describe 
the variation of the amplitude of a d-wave gap, or more 
generally of an anisotropic gap, in direction perpendicu- 
lar to the Fermi surface, and we believe the conclusions 
would be similar concerning the existence of an average 
gap slope. Thus we argue that our analysis is relevant 
irrespective of whether the gap has s-wave or d-wave or 
any other symmetry. 

We note also that recent work |l^] has argued that a 
gap with d-wave symmetry together with an assumption 
of directional tunneling can explain the experimentally 
observed tunneling asymmetries. We will discuss that 
work in the conclusions. 
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In the next section we discuss the tunneling formal- 
ism. Sect. Ill deals with barrier effects, and Sect. IV 
with density of states effects. In Sect. V we consider the 
effect of an energy-dependent gap function. Section VI 
discusses the expected thermoelectric effect in the pres- 
ence of all these effects. We conclude in Sect. VII with 
a discussion. 



II. FORMALISM 

The tunneling probability for an electron of energy E 
above the Fermi level of the tip is given within the WKB 
approximation by 



T(E,V) = e 



-2S(E,V) 



S(E,V) 



E) 



(4a) 



(4b) 



with m the electron mass and V(x) the barrier potential. 
We assume a trapezoidal barrier, with work functions </>t 
and 4>s for tip and sample respectively, and a voltage V 
of the sample relative to the tip, as shown in Fig. I, and 
obtain 



S(E,V) = 



2d 



3(& + V 
- Ef' 2 - ( 







V - Ef' 2 } 



(5) 



for the case where both the arguments raised to the 3/2 
power are positive. When the voltage or energy are such 
that one of the arguments is negative, Eq. (5) still holds 
if that argument is set equal to zero. If the work functions 
<f>t and 4> 8 are similar and the voltage V is small Eq. (5) 
reduces to 



S(E,V) 



V 
~2 



E) 



1/2 



(6) 



with <f> — {(j) s + <f>t)/2 the average work function. The 
derivation assumes V, V + E « <f> and \(f>t — <j) s \ « <t>- 
We assume that electrons tunneling into the supercon- 
ductor go into a band of width D. Within a BCS formal- 
ism the tunneling current from sample to tip is given by 



V is the voltage of the sample relative to the tip. The 
coherence factors are given by the usual form 
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E(e) 



(8a) 



(8b) 



and we assume that the relation between quasiparticle 
energy E and band energy e is given by Eq. (2). p is the 
sample chemical potential. 

The tunneling current Eq. (7) can be written as 



(9) 



I st {V) = / du[f t (u - eV) - 

J — DC 

x g t (u - eV)T(u - eV, V) Ps {uj) 
with p s {uj) the superconducting density of states 

,D/2 

P s{lo)= / de[u 2 (e)S(uj~ E(e)) + v 2 {e)6(LU + E(e))}g s (e) 

J-D/2 

(10) 

Broadening effects [B may be included by replacing the 
(5-functions in Eq. (10) by 



S(uj) -> -- 



7T UJ 2 + r 2 



(11) 



with r a phenomenological broadening parameter. 

We consider for dennitcness a band of width D = 
0.5eV in the sample, a critical temperature T c = 10QK 
and zero temperature gap A(T = 0) = 30meV. The 
temperature dependence of the gap is assumed to follow 
the conventional weak coupling BCS behavior. Figure 2 
shows for reference the tunneling conductance obtained 
from Eq. (7) assuming unit transmission and energy- 
independent normal state densities of states of sample 
and tip, for broadening factors r = and T = 3meV. 
The amplitude of the peaks in dl/dV is proportional to 
u 2 and v 2 for positive and negative voltage respectively. 
Because u 2 (e) = v 2 (2p — e), tunneling characteristics are 
symmetric in the sign of V in the absence of other effects. 



Ist{V) = I ' 2 de[u 2 (e)[f t (E(e) - eV) - /.(25(e))] 

J-D/2 

x g t (E(e) - eV)T(E(e) - eV,V) 

+ v 2 (e)[f s (E(e))~f t (E(e)+eV)} 

x g t (-E(e) - eV)T{-E(e) - eV, V)]g s (e) (7) 

Here, ft and f s denote Fermi functions for tip and sam- 
ple, which may be at different temperatures, and gt and 
g a are normal state densities of states for tip and sample. 



III. BARRIER EFFECTS 

When we allow for energy-dependent transmission as 
given by Eq. (4), the tunneling spectra develop an asym- 
metry that increases as the temperature is lowered, as 
shown in Fig. (3a). By decreasing the work functions 
of tip and sample or increasing the barrier thickness the 
size of the asymmetry is increased. Thus it would easily 
be possible to match the magnitude and temperature- 
dependence of the observed asymmetry. However, the 
sign of the asymmetry thus obtained, larger conductance 
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for positively biased sample, is opposite to what is seen 
experimentally. 

At first sight the existence of this asymmetry may seem 
puzzling, since it occurs for equal work functions of tip 
and sample and for densities of states that are symmetric 
around the Fermi energy for both tip and sample. How- 
ever it is easily understood by inspection of the diagram 
in Fig. 4. Because the gap in the superconducting den- 
sity of states of the sample blocks the transmission in 
that energy range, the transmission for positive sample 
involves electrons that are closer to the top of the barrier 
than that for a negative sample. As the temperature is 
raised and the gap in the density of states closes, this 
asymmetry will disappear. 

Could energy-dependent transmission give rise to an 
asymmetry of the observed sign? It could happen two 
ways. First, if the transmission coefficient for an electron 
were to decrease as its energy increases. While this can- 
not happen in an STM, it could happen in a conventional 
tunnel junction if the Fermi level is in the bandgap of the 
insulating barrier close to the valence band, as discussed 
by Gundlach |l5[ . Second, energy-dependent transmis- 
sion would give rise to an asymmetry of the observed 
sign if the gap structure in the density of states shown 
in Fig. 4 occurred in the tip rather than in the sample. 
Clearly, neither of these scenarios appears to be relevant 
to the case under consideration. 

Next we consider the effect of differences in the work 
function of tip and sample. This can give rise to asym- 
metry in tunneling, and in fact has been proposed to ex- 
plain asymmetries observed in various conventional tun- 
nel junctions If the work function of the tip is lower 
than that of the sample the tunneling conductance in the 
normal state will be larger for a negatively biased sam- 
ple, in agreement with observations. Figure 3(b) shows 
results for a rather extreme case, with work functions 
4> t — O.leV, 4> s = 3eV. While the asymmetry in the nor- 
mal state is of the observed sign, as the temperature is 
lowered the peaks in the conductance exhibit the asym- 
metry of opposite sign, induced by the energy depen- 
dence of the transmission in the presence of the super- 
conducting energy gap. It is found that no combination 
of work functions and barrier thicknesses can give rise to 
an asymmetry of the observed sign in the peaks of the 
tunneling conductance at low temperatures. 

IV. DENSITY OF STATES EFFECTS 

We next consider the effect of non-constant density 
of states and of bandwidth cutoffs. Figure 5 shows the 
resulting conductance when the normal state density of 
states of the sample and/or of the tip has a linear energy 
variation, reflected in the conductance for T — T c in Fig- 
ure 5. This is obtained for an electronic density of states 
in the tip that increases with energy, and/or one in the 
sample that decreases with energy. As the temperature 



is lowered it is seen that the magnitude of the asymmetry 
in the conductance peaks decreases, simply because the 
peaks arise from band energies increasingly closer to the 
Fermi energy. When we include broadening, as in Fig. 
5(b), the asymmetry persists to lower temperatures, but 
there is still a clear progressive reduction in the magni- 
tude of the asymmetry as the temperature decreases. 

Furthermore, when we include the energy dependence 
of the transmission the asymmetry gets further reduced 
and may even change sign, as shown in Fig. 6. In (a) the 
normal state asymmetry, of the sign seen experimentally, 
arises from a sloped density of states of either the tip or 
the sample. In (b), we assume that the band edge in the 
sample occurs very close to the Fermi energy, as shown 
by the dashed line, also giving rise to an asymmetry of 
the sign seen experimentally at high temperatures. |l6| ] 
In the presence of energy-dependent transmission, the 
asymmetry changes sign as the temperature is lowered. 
The work functions assumed here are not particularly 
small, and smaller work functions or larger tip-sample 
distances would make the obtained asymmetry deviate 
even further from what is observed. 

We conclude that asymmetries originating from non- 
constant densities of states and bandwidth cutoffs can- 
not explain the experimentally observed asymmetry. The 
asymmetry originating from these effects invariably de- 
creases rapidly as the temperature is lowered. The asym- 
metry originating in energy dependence of transmission 
discussed in the previous section does persist to low tem- 
peratures but it is of opposite sign to the one observed 
experimentally. Hence we conclude that the experimen- 
tally observed asymmetry must originate in an intrinsic 
asymmetry of the superconducting state itself. 

V. ENERGY-DEPENDENT GAP FUNCTION 

The simplest generalization of BCS theory with con- 
stant energy gap is to assume that the gap depends on 
wavevector k only through the band energy t k . Hence 
the quasiparticle energy is 

E k {e k ) = V(e fc -M) 2 + A 2 (e fe ). (12) 
For energies close to the Fermi energy, 

A(e)=A( M ) + A'( M )(e-M). (13) 

The minimum quasiparticle energy will now occur not at 
eo = A 1 but rather at 

eo = M-A(M)A'0*) (14) 

and is 

A = £(£o) = A(/!) + 0(A' 2 ). (15) 

The BCS coherence factors Eq. (8) are no longer equal 
at the band energy corresponding to the minimum quasi- 
particle energy, eo; rather, 
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U 2 (e ) = -(1-A') 



v 2 (e ) = -(l + A>) 



(16a) 



(16b) 



with A' = A'(/i). For electrons injected into (extracted 
from) the sample, i.e. positive (negative) sample volt- 
age, the tunneling conductance is proportional to u 2 (v 2 ). 
Hence if the slope A' is positive, that is if the gap func- 
tion is an increasing function of electronic energy, the 
tunneling current and conductance will be larger for a 
negatively biased sample, as seen experimentally. We 
define 

^7)+ = height of peak for V samp i e > (17) 

and similarly for V samp i e < 0, and the tunneling asym- 
metry as the difference in the height of the peaks divided 
by the average peak height 



.4 



dl 

dV 



dl 
dV 



r dl 



)+ + #)+]/2 



and we have at low temperatures |]S| 
A(T 0) = 2A'(/x) 



(18) 



(19) 



Figure 7 shows the tunneling conductance for gap slope 
A' = 0.12, with and without broadening. The asymme- 
try is clearly noticeable in both cases and approaches the 
theoretical limit Eq. (19) for the case of no broadening, 
while it remains somewhat smaller for r ^ 0. 

However the effect of the gap slope needs to be bal- 
anced with that of the energy-dependent transmission. 
Expanding the expression Eq. (5) for the transmission 
we obtain 



T(E, V)=T [l + c(E-eV/2)] 



c{eV~ 1 ) = 0.51- 



d(A) 

(evy/ 2 



(20a) 



(20b) 



with <j) the average of the work functions of tip and sam- 
ple. We find for the tunneling asymmetry at low temper- 
atures 



A(T - 0) = 2(A'(/x) 



cA 



(21) 



as T c is approached. For example, for parameters where 
the intrinsic and barrier-induced asymmetry exactly can- 
cel at low temperatures, the intrinsic asymmetry would 
dominate as T approaches T c . Instead, if the gap slope 
scales with the gap itself, its effect on the asymmetry 
will have the same temperature dependence as that of 
energy-dependent transmission. We will assume that this 
is the case here, so that the gap slope approaches zero as 
T — > T c . It should also be noted that this is the behav- 
ior predicted within the theory of hole superconductivity 

& 

Figure 8 shows the combined effects of intrinsic asym- 
metry and energy dependent transmission, for a case 
where intrinsic asymmetry dominates: A' = 0.12, and 
barrier parameters (f> s — <f>t — leV, d — 10A, leading 
to A = 0.044 in Eq. (21). In the presence of also a 
non-constant density of states, the asymmetry seen in 
the normal state persists strongly as the temperature is 
lowered, Fig. 8(b), unlike the cases where no intrinsic 
asymmetry exists (Figs. 5 and 6). 

Eqs. (20) and (21) imply that the observed asymme- 
try should become weaker as the tip-sample distance in- 
creases and the effect of energy-dependent transmission 
increases. There is some evidence for this effect in the 
data shown in Fig. 13 of Ref. 2. The asymmetry for 
the smallest tip-sample distance (largest conductance) is 
approximately 12.5%, and it decreases to about 5—7% 
as the tip-sample distance increases and the current de- 
creases by about a factor of 10. The change in the tip- 
sample distance over that range will be approximately 



Ad(A) 



InlO 



(22) 



From Eqs. (20)-(22) we may deduce the values of the 
average work function and the change in the tip-sample 
distance in the experiment of Ref. 2. Assuming that the 
asymmetry decreased from 12.5% to 6% when the current 
decreased by a factor of 10 and a gap Aq = 30meV yields 



= 0.53eV 



Ad = 3.L4 



(23a) 



(23b) 



which are not unreasonable. The gap slope implied by 
these data is 



A' = 0.062 + 0Mld(A) 



(24) 



For a typical barrier <f) = leV, d = 10A we have c ~ 
beV^ 1 . Hence for a gap Ao = 30meV the tunneling 
asymmetry will be of the observed sign if the gap slope 
is larger than about 0.08. 

The slope of the gap function, A', may or may not de- 
pend on temperature. If A' is independent of T, then its 
effect will become increasingly important relative to the 
effect of energy-dependent transmission as the tempera- 
ture is raised, according to Eq. (21), since A decreases 



where d is the smallest tip-sample distance in the data 
in Fig. 13 of Ref. 2. For example, A' = 0.17 if 
d = 10 A This analysis furthermore implies that the 
observed asymmetry will switch sign if the tip-sample 
distance were to be increased by another 3A, correspond- 
ing to another factor of 10 drop in the conductance. In 
Figure 9 we show the tunneling characteristics for these 
parameters at low temperatures for a range of tip-sample 
distances illustrating this effect. It would be of interest 
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to obtain more detailed experimental information on the 
dependence of the tunneling asymmetry on tip-sample 
separation to compare with the theoretical expectation. 

We conclude from this analysis that it is possible to 
account for the observed tunneling asymmetry in a sim- 
ple way by assuming an energy-dependent gap function 
with a finite slope at the Fermi energy. The existence 
of energy-dependent transmission does not qualitatively 
modify the results, it only implies that the underlying 
gap slope giving rise to the observed asymmetry is larger 
than what one would have inferred assuming the trans- 
mission was energy- independent. To explain the energy 
dependence of the observed spectra above T c it is fur- 
thermore necessary to assume that there is either energy 
dependence in the normal state density of states and/or 
a large difference in the work functions of tip and sample. 
However, these factors would not be able to explain the 
persistence of the asymmetry to low temperatures. Ac- 
curate measurement of the dependence of the tunneling 
asymmetry on tip-sample distance at low temperatures 
should allow for a determination of the effect of energy- 
dependent transmission and hence for accurate determi- 
nation of the gap slope. As we show in the next section, 
further information on the gap slope may be inferred by 
measurement of a thermoelectric effect. 



VI. THERMOELECTRIC EFFECT IN STM 

Further information on the existence of a slope in the 
gap function may be obtained from experiments where 
the tip and the sample are at different temperatures. To 
our knowledge such experiments at cryogenic tempera- 
tures have not been yet been systematically performed 
p8[ , but we are not aware of any fundamental reason 
that would not allow for thermally decoupling of tip and 
sample and keeping them at substantially different tem- 
peratures in vacuum. Assuming there is no heat con- 
duction path between tip and sample, the only way to 
transport heat is through radiation. For example, as- 
sume the sample is at higher temperature than the tip. 
The power transfered from sample to tip per unit area is 
of order 



P = a(T 4 - T t 4 ) 



a = 5.6710 ergsec 



(25a) 



(25b) 



To maintain a temperature difference between tip and 
sample the tip needs to conduct away the heat it ab- 
sorbs sufficiently fast. A typical metallic conductivity is 
k ~ 10' ergsec -1 cmT 1 K^ 1 , so the temperature gradient 
developed in the tip will be 



VT~ -(T 4 -T 4 ) 
fa 



10 



iK 3 



(T 4 - Tt 



(26) 



For example, for a tip of 1cm length at T = at the end 
far from the sample, the end close to the sample would 



heat up to a negligible T t ~ 10~ 3 K if the sample is at 
temperature T s = 100K. A detailed analyisis of heat 
transfer between tip and sample in an STM is given in 
Ref. 19. 

There have been in fact several recent investigations of 
thermoelectric effects in scanning tunneling microscopy 
with normal metal samples P,fL9[ 



22 1 . In these experi- 



ments, performed at room temperature, the tip is heated 
by laser power to a temperature approximately 10K 
higher than the sample. It is noteworthy that in these ex- 
periments the thermopower measured is invariably neg- 
ative, even when the samples are metals with positive 
bulk values for the thermopower. This is presumably 
due to the contribution to the thermopower from energy- 
dependent transmission, as discussed below. 

If we assume that the gap slope A' is positive, as in- 
dicated by the tunneling asymmetry, a positive current 
will flow from the hotter to the colder electrode in our 
case in the absence of other effects, that is, the system 
will have positive thermopower. This is a consequence 
of the fact that the average quasiparticle charge in the 
superconductor, given by 



(27) 



is positive on the average @, unlike the usual BCS case 
with constant energy gap where quasiparticles are charge 
neutral on average. 

For energy-independent transmission and constant 
densities of states the current from sample to tip orig- 
inating from quasiparticles of lowest energy, i.e. with 
E = Ao, is proportional to 

J st = (l-A')[/ t (Ao-eF)-/ s (Ao)] 
+ (l + A0[/ s (Ao)-/ t (A o + eF)] 
= f t (A -eV)- f t (A + eV) 
- A'[/ t (A - eV) + f t (A + eV) - 2/ s (A )] (28) 

and for small temperature gradient this yields for the zero 
current thermoelectric voltage of the sample 



Vo = ^A'^Al 

e T t 



(29) 



hence the thermopower of the junction, Q, is directly 
proportional to the gap slope: 



Ao 



A' 



(30) 



and is positive if the gap slope is positive. Note that the 
thermopower here increases as the temperature is low- 
ered, unlike the usual situation in metals and semicon- 
ductors where it decreases linearly with temperature as 
T-> 0. 

However, contributions to the thermopower will also 
arise from energy-dependent transmission as well as from 
non-constant densities of states. If we consider the cur- 
rent arising only from quasiparticles of minimum energy, 
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Ao, a non-constant density of states of the sample will 
not contribute at low temperatures. Eq. (28) generalizes 
to 

J st = (1 - A')(l + cA )(l + c t Ao)[/ t (A - eV) - f s (A )] 
+ (1 + A')(l - cA )(l - Ct A )[/ s (A ) - / t (A + eV)] 

(31) 

assuming that the voltage V « Ao- Here, Ct is the 
logarithmic derivative of the tip density of states at the 
Fermi energy 

d = g' t (0)/9t(0) (32) 
and we have approximately 

J st = / t (A - eV) - f t (A + eV) - [A' - A (c + *)] 
x [/ t (A - eV) + f t (A + eV) - 2/ s (A )] (33) 

which contributes to the zero voltage thermoelectric cur- 
rent proportionally to 

J st = ^(A' - A (c + c t ))[f t (A ) - /.(A )] (34) 
e 

For the zero current thermoelectric voltage we have for 
small temperature gradient 

V Q = ^[A'-A (c + ct)}^-^ (35) 
e it 

Hence we conclude that the thermoelectric voltage may 
be of either sign depending on the magnitude of the pa- 
rameters c and Ct arising from energy-dependent trans- 
mission and tip density of states respectively. 

In particular, note that the effect of energy-dependent 
transmission is twice as large here as it is in the tunnel- 
ing asymmetry, Eq. (21). This is because the tunneling 
asymmetry occurs at finite voltage, V — A , where the 
effect of energy-dependent transmission is smaller than 
at zero voltage. Hence it is possible to have a tunneling 
asymmetry of the observed sign, reflecting the intrinsic 
gap slope, together with a negative thermopower, domi- 
nated by energy-dependent transmission. 

It is also possible to obtain a negative thermopower 
just from the effect of tip density of states if c t is posi- 
tive, that could dominate over the effect of the intrinsic 
gap slope. However, note that this corresponds to a tip 
density of states that increases with electronic energy, 
and hence it would give rise to tunneling spectra with 
normal state slope that is opposite to the slope that is 
observed experimentally. 

Figures 10 and 11 illustrate some of the effects dis- 
cussed above. Figure 10 shows current versus voltage in 
the presence of a large temperature difference between 
sample and tip, with the tip colder than the sample. 
The thermoelectric effect arising from intrinsic asymme- 
try and from a sloped tip density of states g t is similar, 
as seen in Fig. 10 (curves labeled a and b), and it cor- 
responds to positive thermoelectric power. However the 



dl/dV spectra (Figs. 11a and lib) are very different, 
with the ones corresponding to the sloped gt exhibiting 
an asymmetry that is opposite in sign to that observed 
experimentally and to that originating in the intrinsic 
effect, Fig. 11a. Conversely, a sloped g t that would 
give rise to an asymmetry of the sign seen experimen- 
tally would give rise to a thermopower opposite in sign 
to what is shown in Fig. 10, i.e. negative thermopower. 
Thus experimental results for tunneling spectra together 
with thermoelectric effect could clearly distinguish be- 
tween competing hypothesis of intrinsic asymmetry and 
energy-dependent gt- 

An energy-dependent g s of negative slope would give 
rise to an asymmetry in tunneling of the sign seen exper- 
imentally (Fig. 11c) and to a thermoelectric effect that 
has the same sign as that given by the intrinsic asymme- 
try, but the magnitude of the thermoelectric current is 
much smaller than in the other cases, as seen in Fig. 10 
(curve labeled c). The zero voltage thermoelectric cur- 
rent arising from all these effects is given by 

2 f°° E 

Its = — / dE ; x 

eR J Ao yfW^Aj 

[ ^ + ( c + Ct )E + Cs^^][ft(E)-fs(E)] (36) 

wit c s = g' s (0)/g s (0). The effect of non-constant sam- 
ple density of states is much smaller because the sin- 
gularity at E — A is cancelled for the term involving 
c s . Thus it would be possible to discern the effect of in- 
trinsic asymmetry versus sample density of states on the 
thermopower by calculating the expected magnitude of 
thermoelectric current and comparing with experiment. 

The energy-dependent transmission gives rise to a 
tunneling asymmetry of opposite sign to that seen ex- 
perimentally, as previously discussed, and to negative 
thermopower. However, because its effect on the ther- 
mopower is twice as large as on the tunneling asymmetry 
it is possible to choose parameters so that the tunneling 
asymmetry is of the sign observed experimentally, as de- 
termined by A', and yet the thermopower is negative, 
opposite to what the intrinsic asymmetry predicts. An 
example is shown in Figs. 10 (curve labeled d) and lid. 
However this requires rather large values of both intrinsic 
asymmetry and energy-dependent transmission parame- 
ters. 

Analysis of the thermoelectric effect as function of 
temperature would also allow to clearly distinguish be- 
tween an intrinsic origin and the other effects. Figure 
12 shows the zero voltage thermoelectric current versus 
sample temperature T s for a fixed temperature difference 
between sample and tip, T t — T s = —0.1, with the tip 
colder than the sample. In the presence of only the in- 
trinsic asymmetry (solid line) the thermoelectric current 
would go to zero at T c . Instead, energy-dependent den- 
sity of states of either tip or sample (dashed and dotted 
lines) can give a thermoelectric current of the same sign, 
but it continues to increase as the temperature is raised 
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above T c . The same is true for the current arising from 
energy-dependent transmission, which in addition is of 
opposite sign to that originating in intrinsic asymmetry 
(dash-dotted line). Results for these cases for temper- 
ature gradient of opposite sign, T t — T s = 0.1, are the 
mirror image across the horizontal axis of these results. 

In contrast, the magnitude of the zero current ther- 
moelectric voltage Vo is rather different for both signs 
of the temperature gradient, being larger when the tip 
is colder than the sample. At low temperatures we can 
derive an expression for V that is not restricted to small 
temperature gradients: 

Vo = [A' - A (c + Ct )]Ml [e (ft-MAo - 1] (37) 
e 

which shows that the thermoelectric voltage can become 
very large at low temperatures when the tip is colder 
than the sample. Again, the voltage goes to zero at T c 
if it originates in intrinsic asymmetry and remains finite 
when it originates in the other factors. Results for the 
tip colder and warmer than the sample are illustrated in 
Fig. 13. 

The effect of broadening on the thermoelectric effect is 
shown in Fig. 14. We assume only intrinsic asymmetry 
is present, but the effect is qualitatively the same when 
the thermoelectric effect originates in energy dependent 
transmission or variations in the densities of states. The 
thermoelectric current at low temperatures increases as 
the broadening parameter increases. The thermoelectric 
voltage is reduced in the presence of broadening: the ef- 
fect is most dramatic when the tip is colder than the sam- 
ple, the large thermoelectric voltages obtained in the ab- 
sence of broadening at low temperatures are significantly 
reduced and become comparable to those obtained for 
opposite sign of the temperature gradient. This occurs 
because in the presence of broadening the thermoelectric 
current at low temperatures is dominated by low-energy 
electrons rather than electrons with energies above the 
energy gap. 

Next we consider the situation where both intrinsic 
asymmetry and energy-dependent transmission exist. As 
the tip-sample distance is varied one or the other ef- 
fect could dominate, as previously discussed. Results are 
shown in Fig. 15 and 16. For sufficiently small tip-sample 
distance, assuming tunneling remains ideal, the intrin- 
sic asymmetry effect will dominate at low temperatures. 
Here we have assumed a gap slope A' = 0.24 and an av- 
erage work function cf> — leV. According to Eqs. (35) 
and (30b), the cross-over between intrinsic- asymmetry 
dominated and barrier dominated thermoelectric effect 
at low temperatures should occur for tip-sample distance 
d ~ 15A Indeed, intrinsic asymmetry is seen to domi- 
nate for d = 5A and d = 10A. Note that the cross-over in 
the thermoelectric voltage at low temperatures when the 
tip is colder than the sample (Fig. 16(a)) would be very 
sharp as function of tip-sample distance. Note also that 
in the tunneling conductance instead the cross-over be- 
tween gap-slope dominated asymmetry and barrier dom- 



inated asymemtry would instead occur only for d ~ 30A 
for these parameters, i.e. the intrinsic asymmetry would 
dominate for a much larger range of tip-sample distances. 
Observation of the behavior discussed here would clearly 
evidence the competition between intrinsic and barrier- 
induced effects and allow for the extraction of the intrin- 
sic gap slope. 

Finally we compare in Figure 17 the behavior of tunnel- 
ing asymmetry and thermoelectric voltage as function of 
tip-sample distance, in the presence of intrinsic asymme- 
try, barrier-induced asymmetry, and energy-dependent 
sample density of states. As discussed earlier, intrinsic 
asymmetry dominates in the tunneling asymmetry over a 
substantially larger range of tip-sample distances than in 
the thermoelectric effect. For both cases, broadening in- 
creases the range of tip-sample distances where intrinsic 
asymmetry dominates. 

VII. DISCUSSION 

We have discussed in this paper various factors that 
may give rise to asymmetry in STM tunneling conduc- 
tance experiments. While the various effects discussed 
here can give rise to asymmetries of either sign, we have 
argued that the recent results on tunneling spectra of 
Bi2Sr2CaCu20s+s provide evidence for an asymmetry 
originating in an intrinsic property of the superconduc- 
tor, an energy-dependent gap function. The slope of the 
gap function as function of electronic energy implied by 
the experimental results is positive. We have suggested 
that further data on temperature and tip-sample depen- 
dence of the asymmetry may be able to further support 
this conclusion. 

Furthermore we have proposed that independent evi- 
dence for intrinsic asymmetry originating in a finite gap 
slope would be provided by thermoelectric experiments 
with STM, which should yield positive thermopower in 
certain parameter regimes. Such experiments, always 
yielding negative thermopower, have already been per- 
formed with normal metals, which suggests that they are 
quite feasible. 

The combination of more extensive high quality tun- 
neling spectra and results for the thermoelectric effect 
should be able to determine unambiguously the sign and 
magnitude of the gap slope in Bi2Sr2G'aCu20$+s- The 
gap slope determines the average quasiparticle charge in 
the superconductor and is thus a quantity of fundamental 
interest. The theory of hole superconductivity jll| has 
predicted that the gap slope, and consequently also the 
average quasiparticle charge, is positive for all supercon- 
ductors, and that its magnitude scales with the critical 
temperature. Temperature and carrier concentration de- 
pendence of the gap slope is also predicted by the theory. 

Coffey and coworkers Jr| have recently proposed that 
the observed asymmetry in tunneling conductance dis- 
cussed here is evidence in favor of d-wave symmetry of 
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the order parameter. Their analysis relies on an assump- 
tion of directional tunneling, and the parameters in their 
model are chosen in order to match the observed asym- 
metry. We note that for the case of a half-filled band and 
a band structure with only nearest neighbor coupling, 
that is an electron-hole symmetric system, their analy- 
sis predicts an asymmetry of the observed sign when the 
tunneling direction is parallel or close to one of the prin- 
cipal axis in the plane for both a d-wave and an s-wave 
gap. Even though the asymmetry predicted is larger for 
d-wave, the essential element giving rise to the asymme- 
try in this calculation appears to be the assumption of 
directional tunneling. While this assumption may have 
validity for the case of point contact tunneling it would 
appear not to be applicable to the STM experiments dis- 
cussed here where the tip is mounted perpendicular to 
the Bi-0 layers and tunneling is expected to yield an 
angular average over the ab-plane density of states 0j . 

Recent STM experiments in the presence of a magnetic 
field |H provide further strong evidence for the existence 
of an intrinsic asymmetry associated with the supercon- 
ducting state. The data of Renner et al clearly show 
that when the STM tip is moved from the vortex core to 
a region between vortices the peak for negative sample 
voltages grows faster than that for positive voltages, irre- 
spective of what the sign of the asymmetry in the vortex 
core was. This strongly suggests that the asymmetry is 
directly associated with the superconducting state. As 
we have discussed, the opening of a gap in the absence of 
intrinsic asymmetry would lead to precisely the opposite 
effect, i.e. the positive voltage peak growing faster, due 
to the energy dependence of transmission. Further sup- 
port for this point is provided by the fact that Renner 
et al state that 'the sharpest contrast to map the vortex 
structure on BSCCO is obtained using the conductivity 
at a negative sample voltage V = — A p /e'; (A p = peak 
position. In other words, the peak at negative voltages 
gives the strongest signal on whether the system is nor- 
mal or superconducting. This is related to the prediction 
of the theory of hole superconductivity that photoemis- 
sion should give a much clearer signal of the transition 
to the superconducting state than inverse photoemission 
pij. Similarly one should find a change in the sign of 
the thermoelectric voltage as function of the tip position 
in measurements in the mixed state, with negative ther- 
mopower measured in the vortex core region and positive 
thermopower in regions far from the cores. 

Future work will involve a calculation of tunneling 
spectra as function of position in the vortex lattice within 
a model with finite gap slope such as the model of hole su- 
perconductivity. Possible signatures of a finite gap slope 
in Josephson tunneling are being investigated. It would 
be of great interest to find other experimental signatures 
of a finite gap slope in superconductors. 
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FIG. 1. Diagram of electron energy levels. V is the voltage 
of the sample relative to the tip, E is the electronic energy 
above the Fermi level of the tip. 4>t and 4> a are the work 
functions of tip and sample. 







FIG. 2. Tunneling conductance for a superconductor with 
zero temperature gap A = 30meV and T c — lOOA', for 
T/T c = 0.1,0.25,0.5,0.85 and 1. Constant densities of states 
of sample and tip and energy-independent transmission is as- 
sumed. In the following figures, the same values for gap and 
temperatures are used, (a) F = 0, (b) F = 3meV. 

FIG. 3. Effect of finite barrier height. Barrier thickness is 
d = 20 A In (a), work functions of tip and sample are equal. 
The conductance is symmetric in the voltage above T c , and 
becomes asymmetric as the gap opens. In (b), the conduc- 
tance above T c is asymmetric due to unequal work functions; 
as the gap develops, the asymmetry switches sign. Broaden- 
ing parameter is F = 3meV. 

FIG. 4. Illustration of asymmetry in tunneling resulting 
from energy gap in the sample. For equal work functions, 
the transmission will be higher for positively biased sample 
because the tunneling electrons are closer to the top of the 
barrier. 

FIG. 5. Effect of non-constant density of states. The 
energy variation of the normal state density of states is lin- 
ear, and given by the conductance curve for T = T c . It cor- 
responds to either a sample density of states that decreases 
with electronic energy, or a tip density of states that increases 
with energy, or a combination of the two. The asymmetry de- 
creases as the temperature is lowered, both without (a) and 
with (b) broadening. The thin lines connecting the peaks are 
drawn to illustrate this effect. 

FIG. 6. Effect of sloped density of states (a) and of band- 
width cutoff (b) together with energy-dependent transmis- 
sion. (j> 3 = <j>t = leV, barrier thickness d — 10A. F = 3meV . 
The dashed line in (b) shows the density of states of the sam- 
ple. 

FIG. 7. Effect of intrinsic asymmetry, without (a) and 
with (b) broadening. Gap slope is A' = 0.12. Constant 
densities of states and energy-independent transmission is as- 
sumed. 



FIG. 10. Current versus voltage for sample and tip 
at different temperatures as given in the figure for various 
cases, with constant densities of states, zero gap slope and 
energy-independent transmission except as indicated below. 
Curve labeled a (solid line): gap slope A' = .12; b (dashed 
line): tip density of states of slope mt/gt(0) = — 4eV~ 1 ; 
c (dot-dashed line): sample density of states of slope 
m s /j s (0) = 4eV~ 1 ; d (solid line): barrier parameters 
<t> a = <pt = 0.5eV", d = 20,4, and gap slope A' = 0.36. No 
broadening is assumed. 

FIG. 11. Tunneling characteristics for the parameters of 
the four cases in Figure 10. (a) sloped gap , (b) sloped tip den- 
sity of states, (c) sloped sample density of states, (d) sloped 
gap and energy-dependent transmission. Note that the sloped 
densities of states of sample and tip give rise to tunneling con- 
ductance of opposite slope but to thermoelectric effect of the 
same sign (Fig. 10). Note also that the tunneling asymmetry 
in (d) is of the sign expected from the gap slope, but the sign 
of the thermoelectric effect in Fig. 10 (curve labeled d) is 
opposite because the effect of energy-dependent transmission 
dominates there. 

FIG. 12. Zero voltage thermoelectric current from sample 
to tip for the tip colder than the sample as function of sample 
temperature for fixed temperature difference between tip and 
sample. Constant densities of states, zero gap slope and en- 
ergy-independent transmission is assumed except as indicated 
below. Solid line: gap slope A' = .12; Dashed line: tip den- 
sity of state of slope m t /gt(0) = — 4eV~ 1 ; Dotted line: sample 
density of state of slope m s /g s {0) = 4eV _1 ; Dot-dashed line: 
barrier parameters 4> s = cj> t = leV, d — 7.S4A. F — in all 
cases. 

FIG. 13. Zero current thermoelectric voltage (of sample) 
for the four cases of Fig. 12. The line convention is the same 
as in Fig. 12. (a) Tip colder than the sample. Note that the 
voltage becomes very large at low temperatures. The inset 
shows the behavior of the voltage close to T c . (b) Sample 
colder than tip. Note that the voltages are much smaller 
than in (a). Note also that both in (a) and (b) only in the 
case of intrinsic asymmetry alone (solid line) does the voltage 
go to zero at T c . 



FIG. 8. Effect of intrinsic asymmetry and en- 

ergy-dependent transmission. Gap slope is A' = 0.12, bar- 
rier parameters are (j> s = 4>t = leV, d — 10A. (a) Constant 
density of states, (b) sloped density of states. 

FIG. 9. Tunneling characteristics in the presence of 
non-constant density of states, intrinsic gap slope and en- 
ergy-dependent transmission, at temperature T/T c = 0.1 for 
various tip-sample distances given in the figure (numbers next 
to the curves). The curves are offset increasingly by 0.5 for 
clarity. Broadening parameter is F = 3meV. Gap slope is 
A' = 0.17 and work functions are <f> s = <j> t = 0.53eV\ 



FIG. 14. Effect of broadening on thermoelectric current 
(a) and thermoelectric voltage (b). Gap slope is A' = 0.12, 
constant densities of states and energy-independent transmis- 
sion is assumed. Solid lines are for F — 0, dashed lines for 
r = lmeV and dash-dotted lines for F = 3meV. The ef- 
fect of broadening is to slightly increase the thermoelectric 
current. The thermoelectric voltage is somewhat decreased 
by broadening when the tip is hotter than the sample, and 
substantially so when the tip is colder than the sample. 
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FIG. 15. Zero voltage thermoelectric current (from sam- 
ple to tip) versus sample temperature in the presence of both 
finite gap slope, A' = 0.24, and energy- dependent transmis- 
sion with barrier parameters (j> s — <f> t = leV and tip-sample 
distances given next to the curves. Since the sample is colder 
than the tip, positive thermopower corresponds to current go- 
ing from tip to sample, which occurs when the barrier is thin 
and the intrinsic asymmetry dominates. For temperature gra- 
dient of opposite sign the curves would be the mirror images 
of these across the horizontal axis. F — 0. 

FIG. 16. Zero current thermoelectric voltage (of sample) 
for parameters as in Figure 15 and (a) tip colder than sam- 
ple and (b) tip hotter than sample. Thermoelectric power is 
positive for thin barriers and low temperatures. The inset in 
(a) shows the behavior close to T c . 

FIG. 17. Thermoelectric voltage (of tip) with V = (solid 
lines) and with F — 3meV (dashed lines) at two temperatures 
(labels next to the lines) for Tt — T 3 = — 0.1T C as function of 
tip-sample distance d. Sample density of states is the same 
as in Fig. 5. The data for Y = 3meV are multiplied by 
50. For small d intrinsic asymmetry dominates, for large d 
barrier-induced asymmetry dominates. Note that broadening 
causes the thermoelectric effect to be dominated by intrinsic 
asymmetry over a wider range of tip-sample distances. The 
tunneling asymmetry at temperature T — 0.1T C is also shown 
for T = and T = 3meV (dash-dotted lines labeled A). Here 
again, the asymmetry remains dominated by the intrinsic ef- 
fect for a larger d in the case where F is finite. Gap slope 
is A' = 0.17, work functions are <j> a — <f>t = 0.53eV. Note 
that the range of barrier thickness where intrinsic asymmetry 
dominates is substantially larger for the tunneling conduc- 
tance asymmetry than for the thermoelectric effect. 
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